Thermodynamics of the quantum Ising model in the two-dimensional kagome lattice 
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In the classical Ising model on the kagome lattice, there are macroscopically degenerate classical 
states with exponentially-decayed spin correlation. As the singular quantum perturbation is in- 
troduced by applying the transverse magnetic field, the quantum ground state remains disordered, 
dubbed by " disorder- by-disorder" . Here we compute the temperature dependence of the spin sus- 
ceptibility and the specific heat for difference strength of the field to understand the effect of the 
singular quantum perturbation in this system. The relevance to the ZnCu3(OH)6Cl2 will be also 
commented. 

PACS numbers: 75.30.Gw, 75.40.Cx, 75.40.Mg 



Macroscopically degenerate states often occur in the 
classical level in the anti-ferromagnetic systems that the 
exchange energy can not be perfectly shared by all the 
nearest neighbor spins. It happens naturally if spins are 
placed in the triangular geometry in the two-dimensional 
lattices. In two dimensions, triangles can share the edges 
to form the triangular network or share the corners to 
form the kagome lattice. If we consider the classical Ising 
dynamics in these two systems, both of them show macro- 
scopically degenerate classical ground states. However, 
the spin correlation in these two systems are very differ- 
ent, which determines the fate of the consequent quantum 
ground states when the quantum dynamics is introduced. 
When the transverse magnetic field is introduced by the 
following Hamiltonian, 



(1) 



where = and a'' are the Pauli spin matrices and 
J is taken to be positive, the quantum ground state in the 
triangular lattice favors an spin-ordered ground state, the 
maximally-fl-ippable state, because the spin correlation 
is critical in the classical model. On the other hand, 
the quantum ground state in the kagome case remains 
disordered because the spin correlation is exponentially- 
decayed [l 

In this paper, we focus on the kagome case and com- 
pute its thermodynamic properties with various T. At 
r = 0, the spin susceptibility shows dramatic upturn in 
the low temperature and diverges at the zero tempera- 
ture. For non-zero T, we also see the significant upturn 
but it saturates at the zero temperature. Our calcula- 
tion is the first results on the thermodynamic quanti- 
ties to distinguish the classical disordered state and the 
quantum disorder states. Furthermore, as the macro- 
scopically degenerate classical states lead to the residual 
entropy at T = 0, the quantum dynamics which lifts 
the degeneracy results in additional peak in the specific 
heat in the low temperature at the order of F. Due to 
these properties, our model might be relevant to the re- 
cent discovered ZnCu3(OH)6Cl2, which shows the ab- 



normal upturn in the spin susceptibility and saturates 
at T = 0. There have been several theoretical papers 
trying to describe this unusual property. A straightfor- 
ward interpretation is the presence of the impurity[3j that 
contributes additionally to the susceptibility. However, 
a naive inclusion of the contribution from the impurity 
can explain the data only above 20K[4'|. Different an- 
gles of views by considering the contribution from the 
spin anisotropy and spin-orbital interaction can some- 
what show the upturn but can not show the satura- 
tion at T = due to the "minus-sign" problem in their 
calculation 0, Q. Our calculation advances their results 
because there is no "minus-sign" problem in our model. 
Moreover, our phenomenological approach to understand 
the magnetic systems also opens a new way to explore the 
low-temperature properties that can never be reached by 
the present quantum Monte Carlo technique due to the 
technical difficulty. 

In the following, we use the quantum Monte Carlo 
technique to compute the thermodynamical quantities of 
Eq. (H]) by using the Trotter-Suzuki approximation [3| . To 
calculate the magnetic susceptibility, we have to include 
the perturbed Hamiltonian SH = —hj^i Si — hy J^i Sf ■ 
Discreting the imaginary time direction by n steps, the 
j^th Trotter-Suziki approximant of the partition function 
is given by 

= e-^-[cosh^]^". 
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where i?^ = F^ -f /i^, iV is the number of the lattice sites, 
Ujk are now the classical variables taking only -1-1 and 
-1, and the summation is over the n stacks of the kagome 
lattices. The cluster algorithm is applied along the imag- 
inary time direction. The thermodynamical quantities 
can be obtained by taking the derivatives on Eq. ^ with 



respect to the corresponding thermodynamical variables. 
The results are summarized as the following: 



N 



'X.ZZ 



A-xx 



N 
T 
N 



dp 

+ „cosh2/^ 

n 

dh ' 
,(< (.^^cff^2 



1 aij^ff 2 
> — < — — > 

N dl3 



dhy 



dhl dhy ' 



+ 

T 

N 



hi r2 , PR 

— n an H T tanh 

nri?2 cosh^ M i?3 n 

n 



1 

'^Tcosh^M 

r). 



ar2 ar 



(3) 
(4) 

I 

(5) 

I 

(6) 



where Xzz, Xyy^ ^^'^ Xxx are defined by druz/dh, 
diriy/dhy and dirix/dT respectively, and rax,y,z are the 
magnetization per site. For a fixed temperature, the 
leading term of the error between the Trotter-Suzuki ap- 
proximant of the thermodynamical quantities and those 
obtained from Eq.(IT]) can be shown to be linear in 
1/71. Therefore, our results are the extrapolation of the 
straight line at n — > 00. For each n and an ensemble, 10^ 
Monte Carlo Sweeps are taken. For each set of parame- 
ters, 64 ensemble averages are used. 

Fig.((T]) and ([2]) are the temperature dependence of the 
inverse of the longitudinal susceptibility Xzz per spin for 
various V and different sizes of the system. L in the 
figures is the length of the linear dimension so that the 
number of total spins N = 3L^. The period boundary 
condition is applied in both directions. The temperature 
is in the unit of J and the Boltzman constant is set to be 
1. In high temperature, all curves are fitted well by the 
Curie- Weiss law regardless F and L. As the temperature 
goes down, the Xzz have the abnormal upturn and sat- 
urate at the values depending on the F at T = 0. Spin 
flipping helps reduce the upturn of the Xzz- Since our 
Monte Carlo results for L = 2 is the same as the ones ob- 
tained by the exact diagonalization, we only report the 
exact diagonalization results for L — 2 here. For F = 0, 
the coupling constant in the imaginary-time direction is 
divergent so that we report the result of L = 2 by the 
exact diagonalization only. At finite F, the ground state 
is quantum disorder. L = 4 looks fairly enough for the 
calculation and the size dependence becomes weaker for 
smaller F. Moreover, since the ground state of the F = 
is known as a classical disorder state we believe that 
the result of L = 2 for F = is representative. Fig. ([2]) is 
the magnifying plot around the origin of the Fig.([T]). It 
can be seen clearly that the Xzz of the F = (classical 
disorder) diverges at T = 0, while the ones at the finite 
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FIG. 1: (Color OnIine)The temperature (in the unit of J) 
dependence of the inverse of the P^r spin. L=2, 4, 6 
correspond to N = 12, 48, 108 respectively. The blue straight 
line is the Curie- Weiss fit. QMC is denoted for the result 
from quantum Monte Carlo and ED is the one from exact 
diagonalization 
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FIG. 2: (Color OnIine)The magnifying plot around the origin 
of Fig. IT}. At r = (classical disorder), the longitudinal sus- 
ceptibility diverges at T = 0. At finite F (quantum disorder) 
, they saturate at finite values at T = 0. (Temperature is in 
the unit of J) 



F (quantum disorder) saturate at T 0. This is the first 
result of the magnetic susceptibility thus far that distin- 
guishes the quantum disorder from the classical disorder. 

Experiments measured the magnetic susceptibility in 
the powder samples. In Fig. ([3]) and ([U, we report the 
temperature dependence of the inverse of the powder sus- 
ceptibility X defined by x — l/3(Xa;x + Xyy + Xzz ). One 
apparent feature is that the size scaling becomes even 
weaker. It is probably because x is an averaged quantity. 
Because of the anisotropy, Xxx, Xyyi and Xzz are differ- 
ent from one other. Surprisingly, Xyy is larger than the 
others and the size scaling is small in Xyy On the other 
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FIG. 3: (Color Onlme)The temperature dependence of the 
inverse of the powder susceptibihty x per spin. L=2, 4, 6 
correspond to N = 12, 48, 108 respectively. The blue straight 
line is the Curie- Weiss fit. QMC is denoted for the result 
from quantum Monte Carlo, and ED is the one from exact 
diagonalization. (Temperature is in the unit of J.) 
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FIG. 4: (Color Online)The magnifying plot around the origin 
of Fig.(|3}. At r = (classical disorder), the powder suscep- 
tibility diverges at T = 0. At finite V (quantum disorder), 
they saturate at finite values at T = 0. (Temperature is in 
the unit of J) 



hand, Xyy ^J^d Xxx become equal at F = 0, which is con- 
sistent with expectation. Furthermore, the x shows the 
upturn in the low temperature and saturates at the zero 
temperature for the smaU F cases. In Fig.(|3]), the differ- 
ence between classical disorder and the quantum disorder 
is again clearly seen. The x of the quantum disorder sat- 
urates at T = 0, and the one of the classical disorder 
diverges. 

In Fig. ([5]), we report the temperature dependence of 
the specific heat with various F. At F = 0, a broad peak 
locating at the order of J is clearly observed. In this 
case, there is also a J— function peak at T = represent- 
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FIG. 5: (Color online)The temperature (in the unit of J) 
dependence of specific heat per spin. L=6 corresponds to 
N=108. QMC is denoted for the result from quantum Monte 
Carlo, and ED is the one from exact diagonalization 



ing the macroscopic degeneracy of the classical ground 
states. At finite F, the transverse field lifts the degener- 
acy and therefore results in an additional peak locating 
at the order of F. For larger F, like 0.25 and 0.5, two 
peaks emerges and a spin gap can be clearly seen. For 
F = 0.025 and 0.05, two peaks are well separated. The 
one at the order of J at these two Fs are overlapping 
with the one at F = 0. This indicates that the quan- 
tum perturbation has the dramatic effect on the ground 
states but not the higher-energy states. Unfortunately, 
our numerical resolution is not able to identify whether 
or not it is gapped or gapless for small but non- vanishing 
F. We remark that the states in the degenerate mani- 
fold at F = are connected by the single spin flip, i.e. 
gapless. Whether or not the gaplessness survives with re- 
spect to the quantum perturbation and extends to finite 
F requires more detail study. 

For the application to ZnCu3(OH)6Cl2, there is only 
one parameter in our model to fit with the experiments. 
Although the F term is understood as the external trans- 
verse magnetic field, it can also be regarded as the inter- 
nal quantum spin flipping matrix. In addition, at a first 
glance, considering the Ising dynamics in the spin- 1/2 
system might be absurd for most of the experts. How- 
ever, it has been shown that even if the spin is isotropic in 
the bulk of Fe, Ni, and V, the effective spin interaction be- 
comes anisotropic when the dimension is reduced, namely 
thin film, as long as there is non-vanishing spin-orbital 
interaction 0]. Therefore, that the spin in the spin-1/2 
system is isotropic in 3 dimensions does not guarantee 
in 2 dimensions. In the presence framework, the satura- 
tion seen in the experiments [1, [T^j can be easily under- 
stood. Namely, the quantum spin flip process dominates 
the low-energy physics as long as the Ising dynamics is 
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relevant. 

Our proposal can be tested experimentally by mea- 
suring the magnetic susceptibility on the single crystal 
sample of ZnCu3(OH)gCl2. If the spin flip process is rel- 
evant, experiments should see the shift of the saturation 
as the transverse magnetic field varies. Because the spin 
flip dynamics is enhanced, we expect the saturation shifts 
down as the field increases. In the very recent manuscript 
this effect is observed in the powder sample. The 
applied magnetic field has two components: one is par- 
allel to the hard plane and the other is longitudinal to 
the easy axis. Because the susceptibility is not sensitive 
to the longitudinal field, the decreasing saturation value 
with increasing applied field can be understood as the en- 
hancement of the spin flipping process by the transverse 
field. 

In summary, we compute the temperature dependence 
of the spin susceptibility and the specific heat of the 
quantum Ising model. Different from the Heisenberg 
model in the kagome lattice, the spin susceptibility in 
our model has the intrinsic upturn and saturation at fi- 
nite value. In the specific heat, we observe the two-peak 
structure and the existence of the spin gap for larger F. 
The finite temperature properties in ZnCu3(OII)6Cl2 in- 
dicates that the spin anisotropy and spin-flipping process 
might be dominant in the low temperature. In this case, 
the robust spin dynamics persists to the zero tempera- 
ture might be the reason that ZnCu3(OII)6Cl2 remains 
disordered at the mK range. In addition, there is no 
minus-sign problem in our model even though kagome is 
a highly frustrated system, which advances other meth- 
ods technically to study the low-energy properties of the 



ZnCu3(OH)6Cl2. 
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tion is done by using the facilities in the Supercomputer 
Center in the Institute for Solid State Physics. 



[1 
[2: 

[3: 
[4: 
[s: 
[6: 
[7: 

[8 
[9 

[lo; 
[11 



Electronic address: | chern@issp.u-tokyo. ac.jp] 

R. Moessner, S. L. Sondhi, and P. Chandra, Phys. Rev. 

Lett. 84, 4457 (2000). 

R. Moessner and S. L. Sondhi, Phys. Rev. B 63, 224401 
(2001). 

Y. Ran, M. Hermele, P. A. Lee, and X.-G. Wen, Phys. 
Rev. Lett. 98, 117205 (2007). 

G. Misguich and P. Sindzingre, unpublished, 
arXiv:07041017. 

M. Rigol and R. R. P. Singh, Phys. Rev. Lett. 98, 207204 
(2007). 

M. Rigol and R. R. P. Singh, Phys. Rev. B 76, 184403 
(2007). 

M. Suzuki, Progress of Theoretical Physics 56, 1454 
(1976). 

J. G. Gay and R. Richter, Phys. Rev. Lett. 56, 2728 
(1986). 

T. Imai, E. Nyt ko, B . Bartlett, M. Shores, and D. Nocera, 
unpublished, |cond-ma t/0703141 

O. Ofer, A. Keren, E. Nytko, B. Bartlett, M. Shores, 
D. Nocera, C. Baines, and A. Amato, unpublished, 
cond-mat/0610540 

F. Bert, S. Nakamae, F. Ladieu, D. L'Hote, P. Bonville, 
F. Due, J.-C. Trombe, and P. Mendels, unpublished, 
,.arXiv:0710.0451. 



